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Abstract 

Using full-potential density functional calculations within local density approximation (LDA), 
we predict that mechanically tunable band-gap and quasi-particle-effective-mass are realizable in 
graphene/hexagonal-BN hetero-bilayer (C/h-BN HBL) by application of in-plane homogeneous 
biaxial strain. While providing one of the possible reasons for the experimentally observed 
gap-less pristine-graphene-like electronic properties of C/h-BN HBL, which theoretically has a 
narrow band-gap, we suggest a schematic experiment for verification of our results which may find 
applications in nano-electromechanical systems (NEMS), nano opto-mechanical systems (NOMS) 
and other nano-devices based on C/h-BN HBL. 
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1 Introduction 

Recent experiments [H E] at Columbia Univer- 
sity showing the successful fabrication of hexag- 
onal BN (h-BN) gated graphene field effect tran- 
sistors (BN-GFETs), consisting of both mono- 
layer graphene (MLG) and bilayer graphene 
(BLG) on h-BN gate/substrate which could be 
made arbitrarily thin (down to a mono-layer 
of h-BN), are among the important develop- 
ments in the physics and material science of 
graphene [H [U El [6]. Graphene supported 
on h-BN exhibits superior electrical properties 
with performance levels comparable to those ob- 
served in suspended samples [H El El E]- On 
the other side, in the context of strain engi- 
neering of electronic properties of graphene, an- 
other significant experimental development [8] 
utilizes piezoelectric actuators to apply tunable 
biaxial compressive as well as tensile stresses 
to graphene on h-BN substrate, which allows 
a detailed study of the interplay between the 
graphene geometrical structures and its elec- 
tronic properties. In this scenario, it is useful 
and interesting to explore the effect of biaxial 
strain on the band gap, Fermi- velocity and the 
quasi-particle-effective-mass of a graphene/h-BN 
hetero-bilayer (representing the smallest and 
the simplest form of BN-GFET), which is not 
only feasible but also desirable, since strain- 
engineering [9] of Si, SiGe, Ge has been success- 
fully used in the semiconductor industry to im- 
pressively improve the performance levels of con- 
ventional metal-oxide-semiconductor field-effect 
transistors (MOSFETs), and recently [10] strain 
is seen as a solution for higher carrier mobility 
in nano MOSFETs. 

In this paper, we theoretically explore the ef- 
fect of symmetry preserving homogeneous in- 
plane biaxial strain (which mimics an experi- 



mental condition [8j in which a flexible substrate 
supported graphene can be strained in a control- 
lable way by piezoelectric actuators) on a hetero- 
bilayer (C/h-BN HBL) consisting of an MLG (C) 
and a mono-layer of h-BN, which is the smallest 
and simplest form of BN-GFET. Previous the- 
oretical studies on graphene on h-BN substrate 
[II] and C/h-BN HBL [El [[3] reported the open- 
ing of a small energy band gap in graphene which 
varies with (a) the stacking order, (b) the sepa- 
ration distance of graphene from the h-BN, and 
(c) the externally applied perpendicular electric 
field. Other graphene/h-BN heterostructures 
such as MLG [H] and BLG [E] sandwiched be- 
tween two mono-layers of h-BN have also been 
theoretically shown to have (external) electric- 
fleld-tunable electronic properties. However, a 
detailed study of strain-engineered band gap of 
C/h-BN HBL is not available, although the ef- 
fects of various strain distributions on mono- 
layer h-BN [16J and graphene UM EH ED] 
have been reported recently. Theoretically it 
is shown that the high band gap of mono-layer 
h-BN is strain-tunable [16J, the gap-less nature 
of graphene is robust against small and mod- 
erate deformations [HI EE! EH EQ]. On the 
experimental side, Raman spectroscopy studies 
[El EH [221 E31 El] of graphene reveal that both 
biaxial [8] and uni-axial [ZQ E2 E31 El] strains 
affect the Raman Peaks; the transport proper- 
ties of strained graphene have been investigated 
by depositing samples on stretchable substrates 
[25] . Insights from these studies suggest a de- 
tectable strain-controlled band-gap in a compos- 
ite structure like C/h-BN HBL. This expectation 
has, in fact, been borne out by the results of our 
present study. It is shown that the ground state 
direct band gap (^ 59 meV) of C/h-BN HBL in- 
creases with biaxial tensile strain and decreases 
with compressive strain, whereas the Fermi ve- 
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locity of charge carriers follows a reverse trend. 
Moreover, the variation of computed average ef- 
fective masses (ra*) of charge carriers (electrons 
or holes) with strain is more pronounced with 
increasing trend for tensile strain and decreasing 
trend for compressive strain. These should af- 
fect the carrier mobility, conductivity and optical 
properties of C/h-BN HBL which may be verified 
experimentally for which we propose schemati- 
cally a realistic experimental set-up. 

2 Computational Methods 

Our first-principles study of electronic structure 
of C/h-BN HBL is based on full-potential (lin- 
earized) augmented plane wave plus local orbital 
(FP-(L)APW+lo) method[26], which is a de- 
scendant of FP-LAPW method [27]. We use the 
elk-code [28 j and the Perdew-Zunger variant of 
LDA[29j, the accuracy of which has been success- 
fully tested in our previous work [30J on graphene 
and silicene (silicon analog of graphene). The 
plane wave cut-off of |G + k| max = 8.0 /R m t 
(a.u. -1 ) (Rmt = the smallest muffin-tin radius) 
was chosen for plane wave expansion in the inter- 
stitial region. The Monkhorst-Pack [31 J A:— point 
grid of 20 x 20 x 1 is used for structural opti- 
mization and of 30 x 30 x 1 for electronic cal- 
culations of C/h-BN BLG with simulated ho- 
mogeneous in-plane biaxial strains up to ±15%. 
Our consideration of strains up to 15% is moti- 
vated by recent theoretical calculations [32 — 34] 
as well as experiments [35] which demonstrated 
that graphene can sustain in-plane tensile elastic 
strain in excess of 20%; Kim et al. [25] have mea- 
sured resistances of graphene films transferred 
to pre-strained and unstrained PDMS substrates 
with respect to uni-axial tensile strain ranging 
from 0% to 30%. Although in-plane homoge- 



neous biaxial strains up to ±15% have been 
considered theoretically in [36], the maximum 
compressive strain graphene can sustain is un- 
clear to us. However, for the sake of symme- 
try and as a matter of curiosity, we have con- 
sidered strains up to —15%. The total energy 
was converged within 2/xeV/atom. We simulate 
the 2D-hexagonal structure of C/h-BN BLG as 
a 3D-hexagonal super-cell with a large value of 
c-parameter (c = 40 a.u.) by considering the 
primitive vectors of the unit cell as 

a = (1/2) a e x - (Vs/2)ae y (1) 

b = (1/2) a e x + (Vs/2) a e y (2) 

c ce z (3) 

where a is the in-plane lattice parameter and 
other symbols have their usual meanings. The 
value of a was varied to simulate the application 
of homogeneous in-plane biaxial stress (since |a| 
= |b| = a) to C/h-BN HBL. Figure 1 schemat- 
ically shows the structure of C/h-BN HBL in 
which (a) represents a typical structure, (b)-(d) 
represent top-down views of three different con- 
figurations. 

3 Results and Discussions 

Our calculated LDA lattice constants of 
graphene and h-BN mono-layer respectively 
turned out as ao (C) = 2.445 A and a (BN) 
= 2.488 A, which are in excellent agreement 
with the reported values for graphene (ao(C) 
= 2.445 A in Ref. (11), 2.4426 A in Ref (17), 
2.4431 A in Ref (37)) and mono layer h-BN 
(ao(BN) = 2.488 A in Ref. (16), 2.4870 A in 
Ref (37)) based on plane wave pseudo-potential 
method. Graphene turned out as gap-less, 
while mono-layer h-BN with its valence-band 
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(a) A typical C/h-BN HBL 



(b) Bl 





(d) B3 

Figure 1: Ball-stick-model of C/h-BN HBL with 
top-down views of three different configurations 
Bl, B2 and B3. (a) A typical structure; (b) Bl, 
one C atom (small ball (yellow)) directly above 
the B atom (big ball (red)) and other C atom 
above the center of h-BN hexagon; (c) B2, sim- 
ilar to Bl with B replaced by N (medium ball 
(green)); (d) B3, one C atom is directly above 
the B and the other C atom directly above N. 
(For interpretation of the references to color in 
this figure legend, the reader is referred to the 
web version of this article.) 
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Figure 2: Total energy (E) vs interlayer distance 
d at different a values of C/h-BN HBL(Bl). 
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Figure 3: Variation of E with d at a common 
value of a — 2.47 A for Bl, B2 and B3 configu- 
rations of C/h-BN HBL. 
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maximum and the conduction band minimum 
located at K point of the Brillouin Zone (BZ), 
revealed a direct band gap of 4.606 eV, in 
agreement with previous calculations of 4.613 
eV in Ref. [16J, 4.35 eV in Ref. [37]. Our 
estimated band gap energy of h-BN mono-layer 
is about 23% less than the experimental [38] 
direct band gap energy of 5.971 eV. 

The calculations for the composite structure 
of C/h-BN HBL were done as follows, unlike the 
earlier authors [HI [12] who used the LDA value 
of lattice constant of graphene (ao(C)= 2.445 A) 
for the three configurations of C/h-BN HBL. As 
per the previous studies [HJE], Bl configura- 
tion of C/h-BN HBL (henceforth we call it C/h- 
BN HBL(Bl) or simply Bl) is the most stable 
configuration with equilibrium interlayer spacing 
d (Bl) = 3.22 A at an assumed a (Bl) = 2.445 
A. For Bl, with which we are concerned in the 
present paper, the calculated ground state total 
energies at different d values for five fixed val- 
ues of a ( = 2.445, 2.450, 2.460, 2.470, 2.480 A) 
were obtained as depicted in Figure 2; from this 
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Figure 4: Comparison of energy bands (within 
LDA) of graphene (C) and C/h-BN HBL(Bl) 
near the K point of the BZ. Fermi energy = 
0.0 eV. 
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Figure 5: Variation of Eg of C/h-BN HBL(Bl) 
with 5. Inset shows the variation of Fermi veloc- 
ity vf with 5. 
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Figure 6: Variation of effective mass (ra* 
charge carriers in C/h-BN HBL(Bl) with 8. 
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we estimated the optimized values of a and d as 
ao(Bl) = 2.47 A and do(Bl) = 3.217 A. Simi- 
lar calculations for B2 and B3 are required for a 
comparison of the relative stability of Bl, B2 and 
B3. However, like previous authors [TT1 [T2l [T3] . 
we got the relative stability and minimum in- 
terlayer distance of Bl, B2 and B3 at constant 
a — 2.47 A as depicted in Figure 3, which cor- 
roborate the reported results [TT1 fT2l fT3j . The 
interlayer spacing of B2 and B3 at a common 
value of a = 2.47 A were estimated as g?o(E>2) = 
3.454 A, d (B3) = 3.497 A which agree with the 
reported results [TTI fT2] of d (B2) = 3.4 A and 
d (B3) = 3.5 A. 

The band structures of unstrained graphene 
and unstrained C/h-BN HBL(Bl) near the K 
point of the BZ are depicted in Figure 4 for 
comparisons. As seen in Figure 4, graphene is 
gap-less, whereas C/h-BN HBL(Bl) has a di- 
rect band gap of Eg= 59 meV which is of same 
order of magnitude as that of graphene on a 
hexagonal BN substrate [TT] (Eg = 53 meV) in 
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the same configuration. Therefore, the number 
of h-BN layers below graphene mono-layer has 
little impact on the band structure of C/h-BN 
heterostructure. It is to be noted that Eg is 
under-estimated here because of the LDA. How- 
ever, in experiments [2] involving graphene on 
h-BN substrates, there is no evidence of this 
gap from transport measurements and this ab- 
sence of gap has been attributed to randomly 
stacked graphene on h-BN [2]. The fabrication 
of a graphene/h-BN heterostructure with a de- 
sired alignment seems to be a technological hur- 
dle at present. However, we note that Fan et 
al. [13] have investigated one type of misalign- 
ment of graphene to h-BN layer (along with the 
three aligned structures Bl, B2 and B3 we have 
considered), which they have achieved by trans- 
lating BN mono-layer a distance of \/3a/6 (a is 
the lattice constant of graphene) with respect to 
graphene along the C-C bond orientation from 
a pattern that corresponds to the B3 configura- 
tion of our present study. This misaligned het- 
erostructure with an estimated equilibrium in- 
terlayer distance of 3.4 A is energetically shown 
[13] to be unstable with respect to the most sta- 
ble configuration that corresponds to the Bl con- 
figuration of our present study. Further, this 
misaligned structure is pictorially shown to have 
a small band gap which vanishes when the inter- 
layer distance reaches 4.0 A. It is to be noted that 
the theoretical study in [13] is concerned with the 
effect of interlayer spacings on electronic struc- 
ture of graphene/h-BN hetero-bilayer and there 
is no consideration of biaxial strain effects. We 
also note that since our simulated homogeneous 
biaxial strain preserves the hexagonal symmetry 
of an aligned C/h-BN HBL (say Bl), a misalign- 
ment is not expected to occur in such C/h-BN 
HBL under homogeneous biaxial strain which 
only affect the size but not the shape of the 



hexagonal structure under study. 

The energy dispersion around K point of BZ is 
linear for graphene and therefore the low energy 
quasi-particles mimic the massless Dirac fermion 
behavior [3l @J, i.e., 

E± — Ei? ±hv F k (4) 

where Ep is Fermi energy, vp is Fermi veloc- 
ity of quasi-particles, and (hk) is the momen- 
tum. In case of Bl, low energy dispersion is 
not only linear in k but also parallel to the low 
energy dispersion of graphene, while very low 
energy dispersion within an energy window of 
(-0.06 eV , +0.06 eV) around E F (= 0.0 eV) 
seems quadratic. Since the major portion of the 
energy dispersions around the K point (Dirac 
point) have identical slopes (hence identical vp 
values) for graphene and C/h-BN HBL(Bl), we 
attribute this to one of the reasons for Bl (in 
spite of having a small energy gap) behaving like 
free standing graphene as observed experimen- 
tally. From the linear energy dispersion portion 
of the plot, the calculated average value of vp 
turned out as 0.8 x 10 6 m/s for both graphene 
and C/h-BN HBL(Bl), in agreement with other 
calculated result [13] and in close proximity of 
the experimental value « 10 6 m/s for graphene 
[31 HI [39] . Assuming that the effective mass 
(ra*) of charge carriers of C/h-BN HBL at Dirac 
point given by the relation [13] m* ~ E g /2vp 2 , 
we estimated m* ~ 0.0079m e for Bl (where 
m e is the free electron mass), which is much 
smaller than the effective Dirac fermions mass of 
0.03m e in bilayer graphene [40] . The high value 
of vp (equal to the vp value in graphene), a very 
low value of m* can render C/h-BN HBL pris- 
tine graphene-like electronic properties. How- 
ever, it is to be noted that the observed pris- 
tine graphene-like electronic properties of C/h- 
BN HBL(Bl) can not be attributed to these fac- 
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tors alone. This is due to the fact that exper- 
imentally the graphene/h-BN systems are fab- 
ricated by mechanical ex- foliation and transfer, 
so the stacking between graphene and BN can 
not be pre-determined and relaxed in such tech- 
niques. 

For C/h-BN HBL(Bl), the calculated varia- 
tion of (i) Eg with in-plane homogeneous biaxial 
strain S = (a — ao)/ao is depicted in Figure 5, 
the inset shows the variation of vp with 5; (ii) 
m* with 8 is depicted in Figure 6. 

Our predicted strain-induced modifications of 
Eg, vp and m* should affect the transport and 
optical properties of C/h-BN HBL, which may 
be probed experimentally in a set-up (thought 
of as a specific amalgamation of the reported ex- 
perimental [HEIIE] designs) like the one schemat- 
ically shown in Figure 7. Although further im- 
provements in the design and materials of Fig- 
ure 7 are possible, the important message is that 
such an experiment is feasible with the available 
technology as applied successfully in the experi- 
ments [UEIE]. I n Figure 7, a simple h-BN back- 
gated BN-GFET [H ^ is rigidly fixed on a vari- 
able biaxial strain inducing system (VSIS) like 
the one described in Ref. {£]. By applying a def- 
inite voltage across the VSIS, a biaxial strain of 
definite magnitude can be induced in BN-GFET. 
The current-voltage characteristics of the back- 
gated BN-GFET at different strain conditions 
may be studied to see to what extent the biax- 
ial strain affects the transport properties of BN- 
GFET. Moreover, an optical probe of the direct 
band-gap of C/h-BN HBL may be performed at 
different strained states of BN-GFET in such a 
set-up. This opto-mechanical probe may ascer- 
tain if C/h-BN HBL can be used in designing 
nano opto-mechanical systems (NOMS). 




Au Au 




VSIS: Variable Strain Inducing System 

Figure 7: Schematic view of the proposed exper- 
imental set-up to study the strain-induced mod- 
ifications of the electronic properties of C/h-BN 
HBL. Specific materials and fabrication schemes 
may be found in the experiments [HE]. 

4 Conclusions 

In conclusion, we have studied and compared 
the band structures of graphene and C/h- 
BN HBL (in biaxially strained and unstrained 
states) using density functional theory based FP- 
(L)APW+lo method. From the parallel slopes of 
the energy dispersion bands near K point of the 
BZ of these structures, we estimated identical 
values of Fermi velocity of the low energy charge 
carriers in them. From this and an estimated 
small value of effective mass m* ~ 0.0079m e , we 
infer the pristine graphene-like electronic prop- 
erties of narrow-gapped C/h-BN HBL. However, 
these results do not necessarily reflect the mea- 
sured transport properties of actual devices. In 
particular, the other important factors determin- 
ing the actual device performance include the 
flatness of the substrate, surface phonon energy 
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and scattering, and the trap charges as in the 
oxide substrates. The strain-induced modifica- 
tions of the energy band gap, Fermi velocity and 
effective mass of charge carriers in C /h-BN HBL 
in its most stable configuration was studied in 
detail with output of some testable results for 
which a schematic experimental set-up has been 
proposed. The results, if verified, may find appli- 
cations in future graphene-based NEMS, NOMS 
and other nano-devices. 
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